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ABSTRACT 


An Investigation of thermal per*'urbatlons of Che solar convective zone 
via changes In the mixing length parameter have been carried out, with a 
view toward understanding the possible solar radius and luminosity changes 
that have been cited In the literature. The results show that: (a) a 
single perturbation of a Is probably not the cause of the solar radius 
change give* In ref. 2, and (b) the parameter W = dlnRe/dlnL, can not be 
characterized by a single value, as has been Implied In recent work 
(refs. 2-5) 

INTRODUCTION 


Recent observations of possible changes of the solar luminosity (ref. 1) 
and radius (ref. 2) have spurred theoreticians to try to model the physical 
mechanisms that might produce such changes. One possible mechanism Involves 
thermal perturbations of the solar convection zone. Such perturbations can 
be modeled (to first order) by perturbing the mixing length parameter a(equal 
to the ratio of the mixing length to the pressure scale height), used In the 
standard mixing length theory of convection. In this paper we present the 
results of such an analysis. Other work in this area can be found in 
references 2 to 5. 

METHOD OF ANALYSIS 


The stellar evolution code described by Endal and Sofia (ref. 6) was used 
for this investigation, utilizing the following assumptions and input data: 

(1) Spherically symmetric Sun in hydrostatic equilibrium 

(2) Latest Los Alamos opacities and equations of statk- 
(ref. 7) 

(3) Nuclear reaction rates of Fowler, et al. (ref. b) 

(4) Standard mixing length theory of convection. . 

(5) The use of 700 Interior zones and a small (6xlO”®M<,) 
static envelope, to Insure numerical accuracy. 
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A one solar mass ZAMS model was evolved to the present age of the Sun 
(Td - 4,7x10^ yrs.)* A match to the present luminosity and radius of the Sun 
was found by varying the initial hydrogen abundance, X, and the value of a. 

The adopted values were X = 0.71242 and a * 2.21772. A sequence was then 
calculated such that the time step was slowly decreased until AT = 1 yr. at 
T^. A perturbation of a in the range 0.05% to 4% was then introduced, and the 
full non-linear, time-aepen-lerit evolution followed either for 200 yrs. (with 
AT = 1 yr.) or for 6x10^ yrs. (allowing the time step to increase), at which 
time the normal evolutionary effects dominated the solution. In order to test 
the effect of varying the size of the small time step, similar sequences were 
calculated using timesteps of 0.75, 2, and 5 yrs. No difference in the 
resulting models was seen. The final results are discussed in the next 
section. 

RESULTS 

The results of this analysis are summarized in Table I. We can divide 
the results into three areas, which are discussed below. In the following 
discussion we adopt the following definitions: 


A(x) 


x - x - = change in x. 

new old 


6(x) 


lA(x)*100]/x , , = percent change in x. 
Ola 


(1) Solar Luminosity: As can be seen in Table I, the percentage change 
in Lq is of the same sign as 6(a), and of a much larger value than the 
corresponding 6(R0)(here E solar radius at optical depth 2/3). The 
relation between 6(L^) and 6(a) is linear for 0.75^ 6(a) ^ 4.00, and can be 
written down as: 6(L^) =0.76 * 6(a). The small change in slope at 

6(a) - 0.75 will be discussed later. The characteristic t^me scale to recover 
its Initial value is a thermal time scale of « 6x10^ yrs. 

(2) Solar Radius: In the solar photosphere (x ^ 2/3), 6(r) shows an 

initial change of the same sign as 6 (a), followed by a subsequent relaxation 
toward the new equilibrium radius (larger for 6(a) negative, smaller for 6(a) 
positive) on a much more rapid time scale than 6 (Lq). For 6 (a) » 2, the 
time scale to cross its original value is ^ 200 yrs. For T > 2/3, (r) 

simply shows an immediate relaxation toward the new equilibrium value. As 
Figure 1 shows, 6(Re), and 6(L^) to a lesser extent, changes its behavior with 
6(a) at 6(a) ■ 0.75. The reason for this can be traced to the finite inter- 
polation scheme used for the envelope boundry conditions. Thus the values 
for 6(R^) and 6(L^) are probably not reliable below 6(a) * 0.75. 



(3) W: A parameter used by investigators In the field to characterize 

the changes of radius and luminosity due to any perturbation is 
W = din R^/din • d(Ro)/5(Lo). Figure 2 shows a plot of Wq vs £(a) for the 
data in Table I. Here ■ W fo : the first model after the perturbation is 
Introduced. As noted, the values of W are not reliable for 6(a) <0.75. We 
see that (a) since 6(R.) changes on a much faster time scale than 6(Lg), W 
will change with time (and will change sign, e.g. af ter r;) 200 yrs. for 
6(a) ■■ 2%), and (b) as Figure 2 shows, the value of W. changes with 6(a). 
Hence, there really is no single value for W, as has been implied by earlier 
investigations . 

CONCLUSION 


The main conclusions that can be drawn from this investigation are 
as follows: 

(1) 6(L,) ■ 0.76 6(a), while 6(R«) displays a marked non- 
linear behavior with 6 (a) . 

(2) There is no one value of W which generally characterizes 
thermal perturbations of the solar convective zone. 

(3) Even for 6(a) - 5 , 6(R^) « 2.10"^ while 6(L*) w 3 ! 

Thus for 6(R<j ) ■ 2.5x10'"^ (ref. 2), the extremely large 
6 (Lq) implied shows that a single perturbation of a is 
probably not the cause of the radius change, although 
further work for 6(a) < 0.75, and for a series of random 
a perturbations is definitely indicated. 

(A) Since our results for <5(a) < 0.75 are probably Incorrect, 
a value of W near 0.075 for small 6(a) is not precluded. 

Thus our results show that by taking into account the full 
time-dependent, non-linear behavior of the problem, the 
entire range of W quoted in the literature (refs. 2 to 5) 
may be generated. 
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Table 1 


6a 6R 

o 


.05 

3.818x10"^ 

.1 

7.644x10"^ 

.3 

2.303x10"^ 

.75 

5.802x10"^ 

.8 

5.974x10“^ 

.85 

6.146xjJ ^ 

.9 

6.319x10"^ 

.95 

6.493x10"^ 

1.0 

6.667x10'^ 

1.25 

7.577x10“^ 

1.5 

8.460x10"^ 

2.0 

10.262x10 ^ 

2.5 

12.119x10"^ 

3.0 

1 4 . 033x10""^ 

3.5 

16.038x10"^ 

4.0 

18.047x10“'' 


6R 

W •• — ^ 

6L 6L 

o e 


.0379 

1.007x10 

.0758 

1.008 

.2277 

1.015 

.5701 

1,018 

.6055 

.987 

.6410 

.959 

.6764 

.934 

.7119 

.912 

.7473 

.892 

.9256 

.819 

1.103 

.767 

1.459 

.703 

1.815 

.668 

2.173 

.646 

2.533 

.633 

2.893 

.624 
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Figure 1. The relation between 6(Lj,) and 6(a), 
(solid line) and between 6 (R^) and 6(a), (inter- 
rupted line). The data for 6(a) <0.75 is unre- 
liable (dashed line) . 



Figure 2. The relation between and 6(a). The 
data for 6(a) < 0.75 is unreliable (dashed line). 
The dashed line for a > 3 is only intended as a 
guide for the eye. 



